From adiabatic to non-adiabatic pumping in graphene nanoribbons 
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Non-equilibrium two-parameter pumping transport through graphene ribbons, attached to reser- 
voirs is described. A tight-binding model is solved using Keldysh formalism, and the crossover 
between adiabatic and non-adiabatic regimes is studied. Pumped dc currents through armchair 
ribbons show properties common in two-dimensional systems. The width-dependent dc current in 
zigzag ribbons, reveals that edge states akin to those in two-dimensional topological insulators, do 
not contribute to pumped transport in the adiabatic regime. The interplay between propagating 
and evanescent modes is discussed. 
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The study of transport mechanisms in quantum de- 
vices via dc currents generated by ac fields has become 
an active research area in recent years. Since the theoret- 
ical proposal by Thouless p] -later developed for closed 
and open systems in several works [MZ]-, the application 
of a periodic perturbation to pump dc charge or spin cur- 
rent was achieved in various experimental settings [S1H]. 
Most of these works focused on the adiabatic regime (low 
driving), in configurations with two or more periodically 
changing parameters that yield dc currents proportional 
to the driving frequency 

New insights into pumping have been arising in stud- 
ies on graphene [10 . A two-parameter pumping model in 
the adiabatic regime, based in the Dirac Hamiltonian, ob- 
tained an enhanced pumped current (as compared with 
semiconductor materials), which was attributed to the 
unusual persistence of evanescent modes in the presence 
of Dirac points [11 . Green's function methods, analyzed 
the effects of resonant tunneling in a similar configura- 
tion, and found anomalous behavior in this regime |12j . 
Other works have analyzed pumping in bilayer systems 
[T3] and mono parametric pumping 14J. 

Among the extraordinary properties of graphene, there 
are those arising from confinement effects. It is known 
that finite sample conductances are strongly dependent 
on the boundaries \TE\ HI] . Studies on isolated graphene 
ribbons revealed the existence of unique features in 
the band-structure due to different edge terminations 
[T71 HE] . The most salient aspect is the remarkable width- 
dependence of the energy spectrum with gapped, gapless 
and highly quasi-degenerate bands as the terminations 
are changed. 

On the other hand, in experimental setups, ribbons 
are connected to metallic leads, a fact with significant 
consequences in the context of dc transport 19J. The 
combined effects of finite sizes and contact to leads, pose 
new challenges in our understanding of non-equilibrium 
transport. The aim of the present work is to analyze 



the properties of quantum pumping in graphene ribbons, 
induced by two local ac gate voltages operating with a 
phase-lag, in a setup that includes the effect of metallic 
contacts. The analysis presented below includes a de- 
scription of the adiabatic regime and the crossover to the 
non-adiabatic regime for different edge terminations, in 
a wide range of driving frequencies, voltage amplitudes 
and doping regimes. We show that dc pumped currents 
in armchair nanoribbons (ANR) exhibit properties com- 
mon to other two-dimensional systems. In remarkable 
contrast, zigzag terminated ribbons (ZNR) of appropri- 
ate width, show a vanishing pumped current in the adia- 
batic regime. This is a manifestation of the subtle width- 
dependent band-structure of these systems. These zigzag 
ribbons are akin to two-dimensional topological insula- 
tors, making our results relevant for those materials also. 
Theoretical treatment- We consider the setup of Figure 
[T] with a graphene ribbon along the :r— direction attached 
to metallic reservoirs. The ribbon's dimensions are given 
in terms of its width W — N y (number of rows) and 
length L = N x (number of atoms) in multiples of the in- 
teratomic distance ao = 0.142nm. Two local ac dephased 
harmonic potentials with equal amplitudes drive the rib- 
bon out of equilibrium. The system is described by the 
Hamiltonian: H(t) = H T n,(t) + H res + H tun , where: 

H lih (t) = -t h Y^ ( C U* + h.c.) + 

<r,r'> 

Y, Y^ Eo + VoC0 ^ n ° t + da ^ C rc C rc,, 
a=L,R r a 



(1) 



where th is the hopping amplitude, c\{c r ) is the creation 
(destruction) operator in state r = (x,y), and the sum 
runs over nearest-neighbor sites only. Static barriers of 
height Eq located at sites Yr, and tr separate the rib- 
bon from the reservoirs. Vo, fio and S a correspond to 
the amplitude, frequency and phase lag between the two 
ac potentials applied on top of the barriers. For sim- 
plicity we take 6^=0 and 5r — 5. A tight-binding 
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FIG. 1. (a) Armchair (N y = 3; N x = 8) and (b) zigzag 
(N y = 4; N x — 5) graphene nanoribbons, attached to reser- 
voirs and driven by time-dependent voltages with amplitude 
Vb, frequency f2o and phase difference 5. 



Hamiltonian for a square lattice with constant a' models 
two semi-infinite reservoirs with N' chains and hopping 

t r \ -tires = 2-^a A^kcv ^k n c k c k a • riere K a = (fva,x, ^a 7 y): 

and £k a = — 2t r [cos(k aiX a' ) + cos(k a , y a' )]. The connec- 
tion between ribbon and reservoirs is given by H tun — 
SaSk Q ,r»( w k Q ,r Q cj Ca C I . ci +H.c). In both reservoirs, 
a = R, L, the chemical potential is the same, with fi = 
indicating the neutrality point. 

To obtain the dc current through the ribbon, we use 
the non-equilibrium Keldysh Green's function technique 
[50] as implemented in Refs. [5TJ 155]. This formalism 
produces dc currents as functions of fj,, V and il , at 
and beyond adiabatic and/or perturbative regimes, in- 
cluding details of the band-structure of the ribbon and 
the contacts. The building block is the retarded Green' 
s function 

G* r ,(t,l/) = -iQ(t-l/)({cr(t),ct,(f)}) 

+oo . , 

= y e -ina t / ™ e -iv{t-tr)g rTi ( n}U y (2) 

n— — oo 

The dc current flowing from reservoir a through the rib- 
bon is given by [55] : 

Jdc,u = T 5Z X! / du;Tr{T a (uj + n£lo) x 

P=L,Rn=-oo- J -° c 

G{n,u))Tp{u)& (n,u>)}[fp(fhu>) - f a (Hu + nMlo)],(3) 

where / a (w) is the Fermi distribution function of reser- 
voir a, r a (u>) is the hybridization matrix (obtained af- 
ter integration of the reservoir degrees of freedom) with 
elements T ra . r > a (ui) = 27r^ k w I . cv ,k cl u'k Q ,r' cv <5(w — £k a ). 
All calculations are carried out at zero temperature (fi- 
nite temperature effects will be considered elsewhere) 
and hopping parameters are chosen such that th = t r . 
Furthermore, we set th as the unit of energy, choosing 
Eq = th and explore the range [0 < Vq < E ] for the 
pumping amplitude. We confirmed that maximum dc 
currents are obtained with 5 = tt/2 [H [21] and this 
value for 5 is fixed in the rest of the paper. 

Armchair ribbon - As shown in Figure Hia), these rib- 
bons are connected to reservoirs through zigzag edges. 
Square and honeycomb lattices have N' — N y with 
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FIG. 2. DC current as function of driving frequency fio f° r 
ANR. (a) [(c)] Semi-conducting (N v = 2, iV^ = 20) [metal- 
lic (N y = 4, N x = 20)] ribbon. Solid (violet) and dashed 
(black) lines are for chemical potential fj, — — 1 and re- 
spectively. Vertical dashed (black) lines indicate position of 
resonant frequencies, (b) [(d)] Jd c as function of \i for semi- 
conducting [metallic] armchair. Thin solid (black), dashed 
(red), dotted (blue), and thick solid (green) lines correspond 
to Qo — 0.01,0.3,0.45, and 0.6 respectively. Inset shows cur- 
rent per unit energy vanishing at \x — 0. Other parameters 
are: barrier height Eo = 1, tunneling parameter w c — 1, and 
pumping amplitude Vo = 0.2. All energies are in units of th- 
The phase difference between potentials is set to 5 — n/2. 



equally spaced atoms. Panel (a) [(c)] in Figure [2] show re- 
sults for the pumped current in semiconducting [metallic] 
ribbons, as a function of $7o for two representative values 
of /i. Dashed lines are for fj, = 0, in the gap region for 
semiconducting ribbons [at the Dirac point in metallic 
ones] . As shown, the current vanishes for the frequencies 
explored in the semiconducting ribbon, remaining finite 
in the metallic ones. For larger fi = —1.0th (full lines), 
pronounced maxima and minima reveal the existence of 
discrete electronic levels in the structure. In metallic rib- 
bons, (or semiconducting ones with chemical potential 
away from the gap), we identify the adiabatic regime by 
the linear increase of the current with fio at low frequen- 
cies [4j[7]. For higher frequencies, the current displays 
maxima and minima as function of Q Q , indicating inter- 
ference between different electronic levels [5Tj. Qualita- 
tively similar results were obtained for various values of 
/i 7^ (not shown). 

Panels (b) and (d) show the current versus \x for differ- 
ent values of Qo- In the adiabatic regime (WLq = 0.01th), 
the current vanishes for values of /i within the gap in 
semiconducting ribbons. For metallic ribbons distinct 
peaks appear, indicating the existence of discrete energy 
levels that originate from the linear bands confined by the 
double barrier structure. Values for resonant frequencies 
could be determined from this data by taking the differ- 
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FIG. 3. DC current as a function of pumping frequency fio 
at jj, = for ZNR. (a) Solid (brown): odd (N y — 3) chains; 
dashed (blue): even chains (N y = 4). (b) Jdc as a function 
of fi in the adiabatic regime (ftSlo = 0.01th) for odd chains. 
Solid (brown), dashed (violet), and dotted (green) lines are for 
N y = 3, 5, and 7 respectively, (c) Same as (b) for even chains. 
Dashed (blue), solid (cyan) and dotted (maroon) lines are for 
N y = 4, 6, and 8 respectively, (d) Jdc as a function of n for 
even chains (N y = 4). Thin (black) solid, dashed (red), dotted 
(blue), and thick solid (green) lines for fio = 0.01, 0.2, 0.3, and 
0.4 respectively. Other parameters are same as in Figure [2] 



ence between the positions of any two peaks. The peak 
at fi = in metallic ribbons reflects the current trans- 
ported by all modes up to /i = 0. The inset in panel 
(d) shows that the current per unit energy vanishes at 
the Dirac point because of the Klein paradox [TTJ ITS] , 
At larger and smaller values of /i, charge transport oc- 
curs through many levels with various interference effects 
taking place. Notice that the contact with metallic leads 
plus static barriers produce a particle-hole asymmetric 
current. It is remarkable however, that for a given reso- 
nant frequency, the sign of the current in metallic ribbons 
is unchanged for a wide range of fillings. 

Zigzag ribbon - The connection between ZNR and reser- 
voirs requires to match lattices with different atomic 
spacings. We follow the approach used in Ref. [T5], as 
shown in Figure [lib). The band-structure of zigzag 
ribbons exhibits a characteristic zero-energy band that 
is usually associated with the presence of edge states 
[17] [18] . Several works have revealed width-dependent 
features in the band-structure with a dramatic effect on 
the conductance [T5] [53] . For long ribbons it has been 
shown [TS] that the zigzag/anti-zigzag (or even/odd) ef- 
fect is a manifestation of the different nature of degen- 
erate states at the band center. For even values of N y 
the degeneracy is accidental while for odd values, it is 
protected by the underlying sub-lattice symmetry, being 
quite robust. 

To gain insight into the nature of these peculiar states, 



we show in Figure |3ja) results for the dependence of Jd c 
on the driving frequency fio for odd- (solid line) and even- 
chain (dashed line) ZNR at /x = 0. Similar to metallic 
ANR, the current exhibits alternating sign changes and 
indications of resonant frequencies away from the adia- 
batic regime. Moreover, even-chain ribbons display the 
usual linear dependence characteristic of small ilo- This 
is consistent with the accidental nature of the degeneracy 
of the states at the center of the band, that is removed by 
the confining potential. As a consequence, these ribbons 
exhibit similar behavior to armchair ones. In contrast, 
odd- ZNRs exhibit no current at all at small frequencies, 
consistent with reported results in the resonant tunnel- 
ing regime |12j . This behavior is more clearly exposed 
in the dependence of J^c with /i as shown in panels (b) 
and (c) for ribbons with different widths and fixed aspect 
ratios. Edge states (that persists even in the presence of 
metallic leads) in odd-chain ribbons do not contribute to 
current transport in the adiabatic regime. 

An analysis of the current as a function of [i for dif- 
ferent values of f2o is shown in panel (d), in the case of 
an even-chain ribbon. Solid lines (adiabatic regime) de- 
scribe the characteristic peak structure due to confine- 
ment effects. For resonant and non-resonant frequen- 
cies, the current changes sign more often as compared 
to metallic armchair ribbons, indicating the presence of 
quasi-degenerate states. 

In closing, the expected quadratic dependence of the 
dc current on the pumping amplitude for low amplitudes 
[13 [57], was confirmed for all ANRs and ZNRs studied. 

Role of evanescent modes - Up to here, our results char- 
acterize dc currents due to pumping through propagat- 
ing modes. These are the dominant modes in transport 
when confined systems are long enough. Although con- 
fined graphene shows properties that strongly depend 
on boundary conditions, it is also possible to obtain a 
universal behavior in ribbons under the right conditions. 
The universality is evident when W 3> L and it explains 
the unique values of conductance and enhanced pumped 
currents in the adiabatic limit [TTJ [Jj|] [53] . For shorter 
and wider ribbons, transport is possible because of modes 
from the reservoirs, leaking through the barriers, to the 
system (evanescent). The crossover between these two 
regimes is described in Figure E] Panel (a) shows the 
dc current as a function of length for semi-conducting 
ANR (gap oc 1/W) with the exponential decrease char- 
acteristic of evanescent modes. Panel (b) shows similar 
results for odd-chain ZNR with null current in the adi- 
abatic regime. Thus transport proceeds through evanes- 
cent modes only with an exponentially small contribu- 
tion at long lengths. While the qualitative behavior for 
L <C W and L 3> W is similar to ANR, the the current 
changes sign for intermediate lengths. This crossover pro- 
vides further evidence of the existence of peculiar states 
in these systems. The origin for the sign change in the 
current can be understood as a consequence of the in- 




thus, expected, in the regime satisfying M7q < fc^T. 
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FIG. 4. DC current as a function of length L (=N X ) at /i = 
and £7o = 0.01 for (a) semi-conducting ANR (N v = 14) and 
(b) odd chain ZNR (N y = 19). Other parameters are the 
same as those in Figure [2] 



creasing number of quasi-degenerate edge states as L in- 
creases, with the concomitant opening of new resonant 
frequencies. This means that, while J7 is not a resonant 
for a given L, it can become resonant at a larger one. The 
edge modes that originally did not transport current, be- 
come active in the pumping regime and compete with 
the evanescent modes already present. Although the pic- 
ture presented here captures qualitatively the behavior 
of the current in this intermediate regime, a quantitative 
analysis of the values of the parameters that determine 
the sign of the resulting current, implies further inves- 
tigations of the detailed evolution of the spectrum with 
L. 

Conclusions- We provided a full description of quan- 
tum pumping in graphene nanoribbons with armchair 
and zigzag terminations connected to metallic leads, rele- 
vant for experimental settings. The phenomena observed 
in ANRs is characteristic of non-equilibrium transport in 
2d systems with similar band-structure properties (i.e., 
metallic or semiconducting) [27] ■ In contrast, ZNRs ex- 
hibit a novel width-dependent pumped current that re- 
flects the peculiar nature of edge states. For appropri- 
ate ribbon widths, edge states akin to those existing in 
topological insulator materials have vanishing pumped 
dc current in the adiabatic regime. Furthermore, these 
modes, that persist in the presence of metallic contacts, 
compete with existent evanescent modes originating in 
the contacts. Although the models considered are within 
computationally tractable sizes (few nanometers), they 
correspond to systems available in current experimental 
setups |25| . A simple extrapolation to realistic ribbon 
sizes, predicts dc pumped currents in the range \xA to 
mA at 10 — 100GHz resonant frequencies. Phenomena 
associated with the neutrality point should persist in a 
range of dopings (A/x <C th) as well as at realistic val- 
ues of static barriers and driving voltage amplitudes [5] . 
In general, pumping involves heat dissipation. Within 
the weak driving regime, energy is dissipated at a rate 
ex M7q, while the pumped heat rate is oc UbT^Iq, being 
T the temperature of the reservoirs [5B] . An efficient op- 
eration of the pumping mechanism with low heating is, 



[1] D. J. Thouless, Phys. Rev. B 27, 6083, (1983). 

[2] Q. Niu, Phys. Rev. B 34, 5093 (1986). 

[3] M. Buttiker, H. Thomas, and A. Pretre, Z. Phys. B: 

Cond. Mat. 94, 133 (1994). 
[4] P. W. Brouwer, Phys. Rev. B 58, R10135 (1998). 
[5] B. L. Altshuler and L. I. Glazman, Science 283, 1864 

(1999). 
[6] M. Moskalets and M. Buttiker, Phys. Rev. B 66, 205320 

(2002). 
[7] M. Moskalets and M. Buttiker, Phys. Rev. B 69, 205316 

(2004). 
[8] M. Switkes et. al, Science 283, 1905 (1999); 
[9] P. J. Leek, Phys. Rev. Lett. 95, 256802 (2005); E. R. 

Mucciolo, C. Chamon, and C. M. Marcus, ibid. 89, 

146802 (2002); S. K. Watson et. al., Phys. Rev. Lett. 

91, 258301 (2003); M. D. Blumenthal, et al, Nat. Phys. 

3, 343 (2007). 
[10] K. S. Novoselov et. al, Nature 438, 197-200 (2005). 
[11] E. Prada, P. San- Jose, and H. Schomerus, Phys. Rev. B 

80, 245414 (2009). 
[12] E. Grichuk and E. Manykin, Europhys. Lett. 92, 47010 

(2010). 
[13] C. G. Rocha et. al, Phys. Rev. B 81, 115435 (2010). 
[14] Y. Gu et. al, J. Phys. Cond. Mat. 21 405301 (2009); P. 

San-Jose et. al, Phys. Rev. B 84, 155408 (2011); L. E. F. 

Foa Torres et. al., Appl. Phys. Lett. 99, 092102 (2011). 
[15] D. A. Abanin and L. S. Levitov, Phys. Rev. B 78 035416 

(2008). 
[16] A. R. Akhmerov et al; Phys. Rev. B 77, 205416 (2008). 
[17] K. Nakada et. al, Phys. Rev. B 54, 17954 (1996); L. 

Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 
[18] M. Zarea and N. Sandler, New J. of Phys. 11, 095014 

(2009). 
[19] J. P. Robinson and H. Schomerus, Phys. Rev. B 76, 

115430 (2007). 
[20] L. V. Keldysh, Soviet Phys. JETP 20, 1018 (1965); H. 

Pastawski, Phys. Rev. B 46, 4053 (1992); A. P. Jauho, 

et al , Phys. Rev. B 50, 5528 (1994). 
[21] L. Arrachea, Phys. Rev. B 72, 125349 (2005); %bid 75, 

035319 (2007). 
[22] L. Arrachea and M. Moskalets, Phys. Rev. B 74, 245322 

(2006). 
[23] Z. Li et. al, Phys. Rev. Lett. 100, 206802 (2008). 
[24] J. Tworzydlo et. al., Phys. Rev. Lett. 96, 246802 (2006). 
[25] Z. Chen,' et. al, Physica E 40, 228 (2007); M. Y. Han 

et al; Phys. Rev. Lett. 98, 206805 (2007); L. Jiao et. al, 

Nano Res 3, 387 (2010); N. S. Safron et. al, Adv. Mat. 

24, 1041 (2012). 
[26] L. Arrachea, et al, Phys. Rev. B 75, 245420 (2007). 
[27] T. Kaur, L. Arrachea, and N. Sandler (unpublished). 



